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Abstract 

We consider the statistical properties of the local density of states of a one-dimen- 
sional Dirac equation in the presence of various types of disorder with Gaussian 
white-noise distribution. It is shown how either the replica trick or supersymmetry 
can be used to calculate exactly all the moments of the local density of states. 
Careful attention is paid to how the results change if the local density of states is 
averaged over atomic length scales. For both the replica trick and supersymmetry 
the problem is reduced to finding the ground state of a zero-dimensional Hamiltonian 
which is written solely in terms of a pair of coupled "spins" which are elements of 
u(l, 1). This ground state is explicitly found for the particular case of the Dirac 
equation corresponding to an infinite metallic quantum wire with a single conduction 
channel. The calculated moments of the local density of states agree with those 
found previously by Al'tshuler and Prigodin [Sov. Phys. JETP 68 (1989) 198] using 
a technique based on recursion relations for Feynman diagrams. 
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1 Introduction 



In the presence of disorder mesoscopic systems such as quantum wires and 
quantum dots have electronic properties which are characterised by large sta- 
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tistical fluctuations between different samples [1]. For example, the conduc- 
tance of a typical sample can be quite different from the ensemble averaged 
conductance. A widely studied phenomena, both experimentally and theoret- 
ically, is that of universal conductance fluctuations: the root mean square of 
the conductance is of order e 2 /h (where e is the electronic charge and h is 
Planck's constant) and is independent of the strength of the disorder or the 
sample size [2]. 

Field theoretical techniques based on supersymmetry and the replica trick 
have proven to be extremely useful in studying non-perturbative effects as- 
sociated with disorder in non-interacting electron systems. Supersymmetry 
has been used to evaluate exactly the disorder averaged density of states for a 
one-dimensional Schrodinger equation with a random white-noise potential [3], 
various random one-dimensional Dirac equations [4-7], and a number of mod- 
els associated with the lowest Landau level of a two-dimensional electron gas 
in a high magnetic field [8] . The replica trick has been used to evaluate exactly 
the disorder averaged density of states of the same one-dimensional models 
involving the Schrodinger equation [9], and Dirac equation [10,7]. Non-linear 
sigma models based on either the replica trick [11] or supersymmetry [12-15] 
have also proven to be extremely powerful tools for studying Anderson localisa- 
tion in non-interacting electron systems. Supersymmetric spin chains describe 
the delocalisation transition associated with the quantum Hall effect [16]. 

In a disordered electronic system static fluctuations in the local electric po- 
tential will produce statistical fluctuations in the local density of states. The 
corresponding probability distribution function was first studied by Al'tshuler, 
Kravtsov, and Lerner for the case of two and 2 + e dimensions [17]. Their 
results were based on a perturbative renormalisation group analysis of the 
non-linear sigma model introduced by Wegner [11]. They found that near the 
mean value of the density of states the distribution function was Gaussian 
but slowly decayed with tails of logarithmically normal form. Near the metal- 
insulator transition due to Anderson localisation the distribution approached 
a logarithmically normal form for all values of the local density of states. 
Mirlin and Fyodorov [14,15] have developed a supersymmetric formalism to 
calculate all the moments of the local density of states for a graded non-linear 
sigma model relevant to disordered systems in one, two, and three dimensions 
and a quasi-one-dimensional wire. A novel path-integral approach was used 
by Kolokolov to calculate the probability distribution function for the inverse 
participation ratio of both a finite and an infinite disordered one-dimensional 
system [18]. 

Al'tshuler and Prigodin [19] evaluated exactly all of the moments of the local 
density of states for a one-dimensional metal. They applied the Berezinskii 
diagrammatic technique [20] to the microscopic Hamiltonian. From the mo- 
ments it is straight-forward to construct the complete probability distribution 
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function. They also showed that due to the resulting distribution in Knight 
shifts the distribution function could be related to the line shape of a nuclear 
magnetic resonance (NMR) line. With decreasing temperature and increasing 
disorder the line shape becomes increasingly asymmetric, with large differences 
between the average and typical values. Recently, the results of Reference [19] 
were generalized to the case of a finite-sized ring threaded by a magnetic 
flux [21]. Efetov and Prigodin have used the supersymmetry method to cal- 
culate the NMR line shape for a small metallic particle which corresponds to 
the case of zero dimensions [22,23]. 

Due to advances in nanotechnology it is now possible to experimentally test 
some of the theoretical predictions concerning fluctuations in the local density 
of states. Recent nuclear magnetic resonance experiments on monodisperse 
molecular metal clusters of platinum atoms have tested the theoretical pre- 
dictions for zero-dimensions and find a broad asymmetric line but that the 
density of states fluctuations are obscured by additional contributions from 
surface effects [24]. Experiments have imaged the local density of states in a 
three-dimensional disordered metal (heavily doped GaAs) by measuring the 
differential conductance associated with resonant tunneling through impurity 
states in asymmetric double-barrier heterostructures [25]. Fluctuations in the 
local density of states were detected and found to be enhanced in a magnetic 
field [26]. 

In this paper we show in detail how both the supersymmetry and replica 
trick methods can be used to derive all the moments of local density of states 
for a random one-dimensional Dirac equation. Balents and Fisher [6] recently 
used supersymmetry to evaluate the disorder- average of the one-particle Green 
function for a one-dimensional Dirac equation with a random mass. Bou- 
quet [7] used group theory to simplify their analysis and also treated the 
problem using the replica trick. Our analysis extends their formalism to allow 
for the evaluation of the disorder average of products of Greens functions. This 
is explicitly evaluated for the case of a quantum wire. Our results agree with 
those of Reference [19]. However, in our view the field theoretic method used 
here is more transparent and less cumbersome than the Berezinskii technique. 

In section 2 we introduce the model Hamiltonian, the one-dimensional ran- 
dom Dirac equation and briefly mention the different physical systems it is 
relevant to, ranging from quantum wires to random spin chains. In section 
3 we introduce different forms of the local density of states and define the 
generalized participation ratios that are a measure of the spatial extent of lo- 
calised wave functions. We show how the the moments of the density of states 
and the participation ratios are related to one another. These quantities can 
be expressed in terms of the disorder average of products of retarded and 
advanced Green's functions. Section 4 shows how these products of Green's 
functions may be written in path integral form. Either the replica trick or 
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super symmetry can then be used to evaluate exactly the average over disor- 
der. One is then left with a one-dimensional interacting field theory. Section 
5 shows how a transfer matrix approach can be used to reduce this to study- 
ing a zero-dimensional transfer Hamiltonian. For both the replica trick and 
supersymmetry this Hamiltonian involves a pair of "spin" operators from the 
non-compact algebra su(l, 1). Section 6 shows how the moments of the local 
density of states can be written in terms of the expectation value of compo- 
nents of these operators in the ground state of the transfer Hamiltonian. In 
section 7 we find the ground state for the particular case of the random Dirac 
equation corresponding to a quantum wire. This is then used to evaluate all 
the moments of the density of states. In section 8 we make some corrections 
to the moments of the density of states. In section 9 we use the moments of 
the density of states to calculate the distribution of the density of states. 



2 The model 



The one dimensional random Dirac equation Hamiltonian, h, for a relativistic 
particle with mass M in a magnetic field <3> and in a scalar potential V is a 
2x2 matrix 

h = -i<T z d x + §(x)<T x + M(x)<T y + V(x). (1) 

The Oj's are Pauli spin matrices. For constant <£>, M and V this model describes 
a one dimensional semiconductor when particles are near the Fermi energy 
with a narrow forbidden band. The Hamiltonian is a 2 x 2 matrix because 
it describes particles moving in both directions in one dimension. The wave 
function is of the form \1> = (tp^ip^) where ip^ and ipi correspond to right and 
left moving electrons, respectively. 

The variables $, M and V are random variables with a one dimensional spatial 
dependence. We can separate each of these random variables into an average 
part and a random part 

&(x) = <f> + <j>(x), M(x) — m + rh(x), V(x) — v + v(x). (2) 

The random part of each of the three variables is assumed to have a Gaussian 
white noise distribution centred about zero. If is the disorder strength of 

fa) 

(4>(x)t(x')) = 2D^5{x - x'), (fa)) = 0, (3) 

and similarly for rh(x) and v(x). The disorder strength of rh(x) and v(x) is 
D m and D v respectively. We shall use the general form of the Dirac equation 
shown in equation (1) up until section 7. 
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The random mass Dirac equation in a random magnetic field is when there is 
no scalar potential, V = 0. It also known as an incommensurate Dirac equa- 
tion. The Green's function and localised density of states and has been calcu- 
lated for this model [27,28] . The incommensurate Dirac equation is equivalent 
to the random XX spin chain in a random magnetic field [29,30]. In [6] Ba- 
lents and Fisher consider a random mass Dirac equation where $ = V = 
so only M is non-zero. This is a commensurate Dirac equation [31,4]. The 
commensurate Dirac equation is equivalent to the XY spin chain in a random 
transverse magnetic [29,30]. Bocquet [7] considered a number different cases, 
the random mass model (like [6]), the massless particle in a scalar potential 
(M = $ = 0) and a multiple disorder model (m = = v). 

In section 7 we shall consider the special case of the incommensurate Dirac 
equation. The average parts of M and $ vanish, m = <fi = 0, and there is 
no scalar potential, V = 0. We make the further assumption that the disor- 
der strength of the mass and the magnetic field are equal, = D m . This 
corresponds to a disordered quantum wire with a single channel. 



3 The probability distribution function 



3. 1 The different density of states 



The local density of states at the position x for a general one-dimensional 
system with Hamiltonian h is defined by 

p(E,x) = (x\5(E - h)\x) Q = £ \Mx)\ 2 5(E - ej ), (4) 

3 

with eigenfunction ipj and eigenvalue ej. The subscript Q indicates a quantum 
average (which uses the eigenfunctions as weights) rather than a disorder 
average. For a finite closed sample the electron states make up a discrete 
spectrum. However, the function p can only take on two values, zero or infinity. 
We can define a new and more appropriate density of states from [19] 

/■CO 
p(E',x)f(E-E')dE', (5) 

where / is some arbitrary weight. We define the density of states pi(r),E,x) 
as the density derived from the weight 

/,(*) = -~r—^ (6) 

71 I] 2 + X 1 

where r\ is small, positive and real. This weight is equivalent to the Dirac delta 
function, 5(x), in the limit ?7 — > 0. To obtain the thermodynamic limit of p\, 
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L — > oo, we take the limit rj — > 0. In the thermodynamic limit p(E,x) = 
Pi(0,E,x). 

A full one-dimensional wave function of the general one-dimensional system 
for states near the Fermi energy is related to the spinor ^ = (ip^ipi) in the 
Dirac equation by 

il>(x)=ik(x)e ik " + il> i (x)e- ik ", (7) 

where ip^ and ipi denote the amplitudes for right and left moving electrons, 
respectively. These amplitudes are slowly varying on atomic length scales of 
order kp 1 but the exponential terms vary rapidly over this length scale. The 
intensity of the total wave function is 

mx)\ 2 = \^(x)\ 2 + \^(x)\ 2 + ^(*M(s) e - a *'* + r t Mx)e 2ikFX . (8) 

Thus, the local density of states defined in equation (4) contains terms which 
oscillate rapidly over atomic length scales. We can define an alternative density 
of states that does not include these rapid oscillations 

p(E,x) = J2 \\^(x)\ 2 + l^|(*)fl S(E - e 3 ). (9) 



We can write an equation similar to (5) but with p substituted for p 

/oo 
p(E',x)f(E-E')dE > . (10) 
-oo 

We define pi(i], E, x) as we defined pi(r), E, x), using the weight f v in equation 
(6), but the former is obtained from equation (10) rather than equation (5). 
We shall find, as did Al'tshuler and Prigodin [19], that these two different 
local density of states have quite different probability distribution functions. 

For the Dirac equation the Green's function is a 2 x 2 matrix. A component 
of the retarded and advanced Green's function matrix is defined as 

Ga,p(E,x,x) = (x,a 

where a and (3 are spin indices. In this paper we are interested in the trace 
of this matrix, G R ' A (E,x,x) = Tr spin Gf ]p{E, x, x). We shall refer to G R ' A as 
the retarded or advanced Green's function. In terms of Green's functions the 
density of states is 

Pi(v, E,x) = ^r- (G A (E, x, x) - G R (E, x, x)) . (12) 
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3.2 Participation ratio 



There are a number of different criteria to distinguish energy regions of ex- 
tended and localized states of a particle in a random potential. The inverse 
participation ratio, Pj, is the sum over all positions of the fourth power of the 
jth eigenfunction 

^• = EhM*)l 4 - (is) 

X 

The eigenfunction at position x is ipj(x) with a corresponding eigenvalue tj 
and Hamiltonian h which acts on the wave function \x) as follows 

h\x) = hJ2^j(x) = J2 e A( x )- ( 14 ) 

j j 

The inverse participation ratio is positive for localized states but vanishes for 
extended states in the thermodynamic limit. The participation ratio averaged 
over the ensemble (i.e. the random potential), (Pj), measures those sites which 
make a significant contribution to the eigenfunction normalization. The par- 
ticipation ratio is related to how much the wave functions extend through 
space. For localized states (Pj) is nonzero. If the eigenstates are completely 
delocalized then (Pj) vanishes for an infinite system. 

Instead of considering Pj which is for individual eigenstates of individual sys- 
tems we consider the ensemble averaged quantity in a narrow energy window 
(spectral average) centred about the energy E 

P {2 \E) = (^\^(x)\^(E-e 1 )\(^£5(E-e J )\ \ (15) 

where the angled brackets refer to the average over the ensemble, or the dis- 
order average. Once we take the disorder average of some quantity it will no 
longer depend on the position x. Hence, after summing equation (4) over all 
x and taking the disorder average we obtain 

2L(p(E)) = ^5(E-ej), (16) 

j 

where 2L is the length of the system. Similarly, the average in the numerator 
of equation (15) does not depend on the position, x, so 

P {2 \E)( P {E)) = \^(x)\*8(E - e,)^ . (17) 
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We can generalise equation (17) to 

P^(E)(p(E)) = \^)\ 2k KE - , (18) 

where P^ is the fcth moment of the participation ratio. The first moment 
(k — 1) of the participation ratio is unity. 



3. 3 Moments of the local density of states 



The kth moment of the local density of states, pi, is the disorder average of the 
kth power of the density of states, {pi(rj, E) k ). We can write the kth moment 
of the density of states in terms of the density of states, p, 

/oo 
dE l ...dE k f v (E-E l )...f v (E-E k ) 
-oo 

x(p(E l )p(E 2 )...p(E k )). (19) 

Wegner [11] showed that in the limit rj — > 

( Pl ( V ,E) k ) = (^\^\ 2k S(E - J^dE 1 f v (E-E 1 ) k . (20) 

The disorder average part is the kth moment of the participation ratio defined 
in equation (18) so the kth moment of the density of states is related to the 
fcth moment of the participation ratio by 

/•oo 

(p 1 (r ] ,E) k } = P^(p(E)} dEif v (E - E x ) k . (21) 

J — oo 

We can evaluate the integral over the weight function 

II dxf ^ x)k = Ml ~ fe£ rW^' (22) 

so that 

( Pl ( V ,E) k ) = {^) 1 ~ k n f { ~ )2 l) P [k \ P m- (23) 

From equation (18) P^ = 1, since the left hand side is (p(E)). This agrees 
with the above equation in the thermodynamic limit since in this limit p\ = p. 

The general equation (5) for the density of states pj can be generalised to 
resemble equation (21) [19] 

( Pf (E) k ) = P^(p(E)} / dE 1 f(E - E,) k . (24) 

J — oo 
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The moments of the participation ratio are independent of the weight /. We 
can write a similar equation for the moments of pf by replacing p with p in 
the above equation. The kth moment of p\ can be written in terms of Green's 
functions by using equation (12) 

(P^,E) k ) = j^-y- g J0^y_ (G A (E, x, x) k -iG R (E, x, xy) . (25) 

In this paper we calculate the Green's function correlator 

(G A (E,x,x) k - j G R (E,x,x) j ) . (26) 

This correlator enable us to calculate the moments of the density of states, 
(p\). By extending this method a little we can also calculate (p k ). This allows 
us to calculate the moments of the participation ratio via equation (23). 

In section 7 we will calculate the moments of the local density of states, 
(pi(r), E) k ), for a random mass Dirac equation. In section 8 we find (pi(r), E) k ) 
for the same model which gives the moments of the participation ratio, p( k \ 
Knowing the moments of the participation ratio we can solve equation (24) 
for (p k ) using a general weight, /. In section 9 we take / to be the derivative 
of the Fermi-Dirac distribution function. 



4 Path integral form of products of Green's functions 

A component of the advanced and retarded Green's function matrix at some 
position x and energy E is defined in equation (11). Note that the retarded 
Green's function is the complex conjugate of the advanced Green's function. A 
solution for the Green's function can often be obtained by writing it in terms 
of a path integral. There are two standard ways to construct the path integral, 
the replica trick and supersymmetry. Both techniques allow the average over 
the disorder to be performed analytically. 



4-1 Replica trick 

To construct the replica trick path integral form of the advanced Green's 
function we firstly write it as a path integral over two complex variable, s and 
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G A (E,x,x) = i 



J Vs*Vss*(x)s(x) 



X 



x exp y—i J dys*(y)[E — irj — h]s(y) 
J T>s*T>sT> exp ( — % J dys*(y)[E — irj — h]s(y) 



(27) 



The region of integration (-L, L) is over the length of the system. We shall 
later take the limit L — > oo. The retarded Green's function is the complex 
conjugate of the trace of the advanced Green's function 



G R (E, x, x) = -i 



J Vs*Vss*(x)s(x) 



X 



x exp \ i J dys*(y) [E + irj — h]s(y) 
J Vs*VsVexp (i J L L dys*(y)[E + ir] - h}s(y) 



(28) 



The dominator is Z, the partition function. 



If the variables s and s* are replicated n times the advanced Green's function 
becomes 



G (E,x,x) = - 
n 



/n 
1=1 

x exp i^-i J L d vYl s ii E - if] - h]sij 
: Jvs*Vsexp^—i J dy ^ s* [E — irj — h] si 



■ (29) 



The dominator is now Z n . Thus, if we let n — > the denominator is unity and 



G A (E,x,x) = lim- 



/n 
Vs*VsY,s:{x) Sl {x) 

x exp (-i dy [E-iri- h]s^j . (30) 



This is one form of the replica trick path integral. A similar expression may 
be constructed for the retarded Green's function. There are other forms of the 
replica trick path integral which only differ in the placement of the complex 
number i, such as the form used by Bocquet [7]. 
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We shall now construct a replica trick form of a multiple of Green's functions. 
The jth power of the advanced Green's function is 



G A (E,x,x) j 



x 



J Vs*Vs (s*(x)s(x)y exp ^— i J dys*(y)[E — irj — h]s(y) 
J Vs*Vs exp ( — % J dys*(y)[E — irj — h]s(y) 



, (31) 



and the jth multiple of the retarded Green's function is the complex conjugate 
of this. The product of the (k — j)th power of the advanced Green's function 
and the jth power of the retarded Green's function may be written as a path 
integral over four complex variables, s, s*, s' and s'*, 



G A (E,x,x) k - j G R (E,x,x 



x 



(k-j)\j\ 

J Vs*VsVs / *Vs / (s*(x)s(x)) k ~ ] (s'*(x)s'(x)y 
exp (-i dy [s*[E -irj- H]s - s'*[E + irj- h]s'^J 

Vs*VsVs'*Vs' 



x 



exp ^-i J L dy [s*[E - ir, - h]s - s'*[E + ir] - h]s']j 



• (32) 



If we replicate s n times and s' n' times then 



G A {E,x,x) k - j G R {E,x,xy = 



i k (-iyr(n)r(n') 



X 



T(k-j+n)T(j+n>) 

r / n \ k ~j In' ^ J 

J Vs*VsVs'*Vs' r£8;(x) 8l (x)\ \'Zs?(x)s' l (x) 



,1=1 



x exp I — i J dy 



J2 s*[E - ir, - h] 3l - £ tf[E + iri- h]^ 
i=i i=i 



x exp — i 



1 l-L ^ 



J Vs*VsVs'*Vs' 

n n' 

s:\E- in- h] Sl - ]T sT[E + h]s\ 



1=1 



1=1 



• (33) 



The prefactor T(n)T(n')/T(k— j + n)T(j + n') is derived in appendix A. When 
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taking the limit n, n' — > 0, as before, the denominator is unity and 



G A (E, x, x) k ~ j G R (E, x, x) j = lim 



i k (-iyr(n)r(n') 

T{k-j + n)T(J + n') 

k-j 



x / Vs*VsVs'*Vs' (j2s*(x) Sl (x)j lj2s' l *(x)s' l (x) 



< exp I — / / dy 



J2 s*[E -ir,- h]s t - J2 s' t *[E + irj- h]s[ 
1=1 1=1 



• (34) 



4-2 Supersymmetry 



The supersymmetric path integral form of the advanced Green's function is 
an integral over two complex variables, s and s*, and two Grassman variables, 
f and T [32,12], 

G A (E,x,x) = i j V£*V£Vs*Vss*(x)s(x) 

x exp (-i J L ^ dy [C[E -irj- h]£ + s*[E - irj - h]s]j . (35) 

Unlike the replica trick no denominator corresponding to a normalization (or 
partition function) appears because the unwanted contributions from the com- 
plex and Grassman variables cancel. The Grassman variables are introduced 
for this reason. As stated before, the retarded Green's function is the complex 
conjugate of the advanced Green's function. The jth multiple of the advanced 
Green's function is 



G A (E, x, x) j = ^ f VCV£Ds*Vs{s*{x)s{x)y 

exp ^-i J^dy [C[E -irj- h]£ + s*[E - irj - h]s^j . (36) 



3- 

x 



The product G A (E,x,x) k ^G R (E,x,xy can be written in terms of a super- 
symmetric path integral over eight variables [12] 



G A {E,x,x) k - j G R (E,x,x) j 



i k {-iy 



(k-j)\j\ 

x J VCVZVs*VsV£*VZ'Vs'*Vs'(s*(x)s(x)) k - j (s'*{x)s'(x)y 

x exp ^-i j dy [l~*[E - irj - h]£ + s*[E - irj - h]s 

-e*[E + irj- h]? - s'*[E + h]s']) . (37) 
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5 The transfer Hamiltonian 



We now briefly review the transfer matrix method [33,3-5] which can be used 
to reduce a one-dimensional field theory into a zero- dimensional Schrodinger- 
type equation. We wish to calculate the disorder average of a product of 
Green's functions, (G A (E, x, x) k ~ : 'G R (E, x, xY) . We have shown in the pre- 
vious section that this average is of the form 

(<f ) (x)}s = Q'vX ( f>(x)e- s my (38) 

where S is some action and a; is a spatial variable. For the replica trick 

and for supersymmetry 

<f>(x) = (s*(x)s(x)) k - j (s'*(x)s'(x)) j . (40) 
The action is defined as the integral of the Lagrangian 

S = [ L Cdy. (41) 



All the disorder is contained in the Hamiltonian, h, which is contained in the 
exponential term e~ s< - x \ If we have a random Gaussian variable, V, which 
has an average value v and disorder strength D v and g is some non-random 
function then the average may be taken by using [4,5] 

/exp (i dyV(y)g(y)\ \ = exp (-D v dyg(y) 2 ) . (42) 



Once we have taken the disorder average we can find a Hamiltonian, known 
as the transfer Hamiltonian [3-5], from the disorder averaged Lagrangian. If 
we have a coordinate q and a momentum p the transfer Hamiltonian is 

H = C-pq. (43) 

A path integral is a sum over all possible configurations. We can rewrite the 
path integral of ((f)) s as such a sum. The continuous interval (—L,L) can be 
written as a regular lattice with spacing a so that there are 2L/a + 1 sites. 
The value of the variable X at site x is X x . Because of periodic boundary 
conditions X_l = Xl. The path integral is now an integral over the variables 
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X x with x = —L, — L + a, . . . ,L. Or more compactly 

(<j)(x)}s = J dX L dX x K(-L,x,X L ,X x )cf>K(x,L,X x ,X L ), (44) 

where most of the integrals are hidden inside the functions K which may be 
thought of as propagators. Note that the inside the integral longer has any 
spatial dependence because we have already taken the disorder average. 

From a propagator such as K we can define a set of eigenfunctions, ^ m L and 
ty m R. The two eigenfunctions, ^ m L and ^ mR , are left and right eigenfunctions 
respectively. We need to distinguish between the two because the transfer 
Hamiltonian need not be Hermitian. The right eigenfunctions and propagator 
should satisfy an equation of the form 

ttmflPQ = / dX y K{y,x,X y ,X x )* mR {X y ). (45) 

Using Feynman's transfer matrix technique [33] we take x = y + a and find a 
Schrodinger-type equation for ^ mR : 

%d ^f = H * mR > (46) 

where H is the transfer Hamiltonian defined in equation (43). The left eigen- 
function satisfies a similar equation. We can define some eigenstates, E m , such 
that H^f mR = E n ^ mR and ^ m LH = E n ^ mL . We define the eigenfunctions to 
be orthonormal so that 

J dXV* mL (X)V nR (X) = 5 mn . (47) 



A propagator can be written in terms of its eigenfunctions and eigenstates 

K(y,x,X y ,X x ) = ^rn R (X y )** mL (X x )e W [-E m (x-y)}. (48) 

m 

We substitute this into equation (44) and using the orthonormality of the 
eigenfunctions obtain 

(0O*O)s = E / \^m(X)\ 2 <pdXexp[-2LE m ], (49) 
where we have defined \^ m \ 2 = ^* m L^mR- 

We shall now derive an equation of the above form for the Green's func- 
tion correlator (G A (E,x,x) k ~ : 'G R (E,x,xy). We obtain both a replica and a 
super symmetric version. The equation can be further simplified in the limit 
L — > oo. 
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5.1 Replica trick 



Since the Hamiltonian, h, for the random Dirac equation (1) is a 2 x 2 matrix 
the path integral variables in equation (34) are two component vectors. For 
example 



(50) 



By substituting the Dirac Hamiltonian (1) into equation (34) we can obtain 
the action 



-/_> 



E [s*i<? z dsi + i^s*a x s t + iMs^a y s + (iV - e_)s*sj] 



.1=1 

n' 



- ]T [sfads't + i&sfasl + iMs' l *a y s' l + (iV - e+)s|*s|] 
i=i 

The coordinates are 

qi = si = (s ih s u ), q{ = s[ = (s' ip s' n ), 
and so we can calculate the momenta 

Pi = (sih s n), Pj = (SiT> s u)- 



(51) 



(52) 



(53) 



To obtain the transfer Hamiltonian the action needs to be transformed into 
the form of a coherent path integral [6,7] 



5 «T 



bii, 



(54) 



and bi = (bi^b^). The transformation of the primed terms is identical. The 
coordinates and momenta in terms of these new variables are 



9i = (&JT> 6 ii)> Pi = (-&iV 6 u)> 
Qi = ( 6 it> 6 J)> Pi = (6?T,-^)- 



The transformed action is 



5 = 



+M(6f T &; i - 6 U 6 IT ) + {iV - e_)(&f t 6, T + b n bl) 

n 



(55) 



+M(6; t T 6; t i -6; i 6; T ) + ^- e+ )(6; T T 6; T + ^;|) . (56) 



.tj.' 
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At this stage we take the ensemble average of the action using equation (42). 
After averaging over the disorder we can find the transfer Hamiltonian [3-5]. 
From equation (43) the transfer Hamiltonian is 

H = -2icf)(Z y - Z' y ) + 2im(Z x - Z' x ) - 2i{v - E){Z Z - Z' z ) + 2 V (Z Z + Z' z ) 
+ AD^Zy - Z' y ) 2 + AD m (Z x - Z' x ) 2 + AD V (Z Z - Z' z ) 2 , (57) 

where 

n n 

Z = -U2pi(ia x ,ia y ,a z )qi, Z' = \ J2p'i(ia x , ia y , a z )q[ 
i=i i=i 

n n' 



i=i i=i 

Zy = \ f>Wl + Z' y = \ ^(bX + b' n b' n ), 



n ii 

\ E(^t + Wl), z' z = \ + b 'n b 'l)- (58) 



1=1 



Note that the transfer Hamiltonian does not depend on the spatial coordinate 
x making it a zero-dimensional equation. We can define some ladder operators 

Z± = Z y t iZ x , (59) 

and similarly for Z'±. The vector Z satisfies su(l, 1) commutation rules [7] 

[Z+, Z_] = -2Z Z , [Z z , Z + ] = Z + , [Z z , Z.\ = -Z-, (60) 

and also (Z' + )^ = Z'_. Identical rules govern Z. 

We define the mth eigenstate of the transfer Hamiltonian as E m with corre- 
sponding eigenfunction \l/ m . Since the transfer Hamiltonian need not be Her- 
mitian the left and right eigenstates may be different. After undergoing the 
transformation in (54) the disorder average of equation (34) in terms of E m 
and ^> m is [3-5] 

G A (E, x, x) k ~ j G R (E, x, xy) = hm ^M)^ W r K) 



n,n'^o T(k -j + n)F(j + ri) 

(n \ k~3 ( n' \ ^ 

£(6fr6iT + biMM rZ(bX + b'nb'l)] 

x-£\*m(b,b\b',b*)\ 2 e- 2LE ™. (61) 

m 

The partition function for n, n' — > is defined as 

Z = lip f Vs*VsVs'*Vs'e- s , (62) 
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and when deriving equation (34) we showed that Zq = 1. After transforming 
as shown in equation (54) the disorder averaged partition function in terms of 
E m and ^ m is 

(Z }= lim [db^Mb*dl/Y l \iS! rn (b,tf,b',lfl)\ 2 e- 2LEm . (63) 



As the system becomes infinite in length, that is L — > oo, the exponential term, 
e -2LE m w jjj vanish for all positive energy levels. In the the disorder averaged 
partition function the dominant term is the ground state energy term 

(Z )= lim / dtfdbdb*dl/\Vo(b,tf,b',tf)\ 2 e- 2LEo (64) 

n,n'— »0 J 

which may alternatively be written as 

(Z )= hp L ^\e~ 2LE ^) R (65) 

n,n — »0 

where and \^) R are the left and right ground state wave functions re- 
spectively with ground state energy E . As stated above this average partition 
function must be unity so E = and 

lim L {*\*) R = 1. (66) 

n,n — >0 

Hence, to calculate the Green's function correlator in the limit of an infinite 
system we need only calculate the ground state of the transfer Hamiltonian, 
H . After taking the L — > oo limit of (61) and using equation (58) to write it 
in terms of Z and Z' 

(G A (E x x) k ~ j G R (E x xY) - lim (2Q*(-lffln)r(n') 

x L (*\Z k -iZ'i\$) R . (67) 

From equation (25) the kih. moment of the local density of states, p±, is 

E) > ~ „teo ^ g jl(k-j)W(k-j + nnj + n>y (68) 
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5.2 Supersymmetry 



We substitute the Dirac Hamiltonian (1) into equation (37) to get the super- 
symmetric action 

S = ~f L dy [?a z dt + s*a z ds - £*a g d? - s'*a z ds' 

+ i$(£V*£ + s*a x s - £'*(T X Z' - s'*a x s') + iM(£V„£ + s*a y s - - s'*a y s') 

+i(V -{E- i v ))(?S + s*s) - i(V -(E + iri))(t? + ■ (69) 

If we define the coordinates as 



the momenta are 



q 1 = s = (s< [ ,s l ), 52 = f = (£t>£i)> 

^ = S ' = ( S ' T ,4), ^ = = (70) 



pi = K,-«i), ^=(^,-c)- ( 7i ) 



As in the replica case we need to transform this action. We shall use the same 
transformation as Bocquet [7] rather than the transformation used by Balents 
and Fisher [6] 



s T — > £ T — > /j, 

*i ^ *!> - /|. 

s-f — > 6|, £j — > /p 

*; - & i> c - -/i- (72) 



The negative sign in the transformation of the Grassman variables is because 
for any Grassman variable (£*)* = — £. In terms of the new variables the 
coordinates and momenta are 



gi = (b h b\), 
Pi = (-^{' & l)' 



ft = (/T./J). 

92 = (/f,tf), 
P2 = 

P 2 = (/?,/{)■ (73) 
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The transformed action is 

S = - £. dy [fjdfr + fM} - ffdfi - f{dfi 

+ &{<% T - bflbl - tfdb\ + b\db'l 

+ mUl + fi/i - fUi - f\f{ + b \ b \ + hh - bpi - b\b\) 

+ M{f\f\ - / t / A - fffi + f\f{ + - b^ - + b\b\) 
+ «(V - - ^))(/ T Vt + /J/i + &J&t + %) 

-i(V ~(E + ir,))(tff{ + fU'l + + • ( 74 ) 



We take the average over the disorder as described in equation (42). Following 
Bocquet [7] we define two superspins with three components 

J = -\v\{io x , ia y , a z )q x - \p 2 {ia x ia y , a z )q 2 , 
J' = \ V \{io\io\o z )q\ + lp' 2 (ia x ,ia\a z )q f 2 . (75) 

Both these superspins satisfy su(l, 1) algebra like Z' and Z in equations (60). 
We can write the transfer Hamiltonian in terms of components of these two 
superspins 

H = -2t<p(J y - Jy) + 2im(J x - J' x ) - 2i{v - E)(J Z - J' z ) + 2 V (J Z + J' g ) 
+ AD^Jy - Jy) 2 + AD m {J x - X) 2 + AD V (J Z - J' z f. (76) 

Note that the transfer Hamiltonian, like the transfer Hamiltonian of the replica 
trick, is independent of the spatial coordinate x. We can decompose the su- 
perspins, J = J + S, into fermionic, S, and bosonic, J, parts 

J = ~^Pi(i(r x ,i(T y ,a z )q u 3' = \p^{ia x , ia y , a z )q[, 

S = -\p 2 {io x , ia y , a z )q 2 , S' = \p' 2 {ia x , ia y , a z )q' 2 . (77) 

The J, J', S and S' all satisfy the su(l, 1) algebra. In terms of b, b', f and /' 



Jx 




J' x = 


WM-b\b\), 


Jy = 


\(b\b\ + b&), 


J'y = 


\{bpl + b\b',), 


J Z = 




J' z = 


\{bp, + b\b't), 


s x = 




$x = 




Sy = 


!(///! + /i/ T ), 


Sy = 




s z = 


|(/t/t + /i/i)> 


s' z = 


UffX + fif?)- 



We can construct an integral like equation (61). The supersymmetric version 
of this equation has, in addition to the integrals over the complex or bosonic 
variables, b, b', and b'\ integrals over the Grassman or fermionic variables, 
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/, /', p and f^. The mth eigenstate of the transfer Hamiltonian, \I/ m , has a 
dependence on the fermionic variables as well as the bosonic variables. Like 
the replica trick we can define a disorder averaged partition function which is 
set to unity. When we take the limit L — > oo we obtain the supersymmetric 
equivalent of equation (67) 

(G A (E,x,x) k ~>G R (E,x,xy} = { ^ k }~^ L Wt j J'JmR, (79) 

where, like the replica trick, \^)r and l are the right and left ground state 
wave functions of the transfer Hamiltonian with ground state energy Eq = 0. 
From equation (25) the kth moment of the local density of states, pi, is 



The first moment of this local density of states is 

(p 1 (ri,E)) = ± L {*\J z + J , z \*) R . (81) 

71 

If we take the limit 77 — > this equation should be equivalent to equation (9). 
Using equation (78) we can relate ip-^ and ipi to the bosonic variables b and b' . 
To find the moments of the full local density of state, pi, rather than pi, we 
must include the rapidly oscillating terms in (8) . This means that instead of 
J z in equation (80) we have J z + \ J + e~ 2tkFX + ^J_e 2lkpX , and similarly for J' z . 
So, 

b\ k 1 

x (J' z + \j' + e- 2ikFX + \.J'_e 2ikFx y\<S>) R . (82) 

Note that when this is expanded out terms involving J + J'_ and J^J' + can 
produce contributions to the moment that are slowly varying in space. This 
argument to find (p k ) is also valid for the replica trick. We used supersymmetry 
rather than the replica trick because the notation is simpler. 



6 The ground state 



6.1 Replica trick 



Bocquet [7] has calculated the ground state for the transfer Hamiltonian as- 
sociated with the one-particle retarded Green's function. This ground state is 
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a linear combination of the basis states {\p}} 

oo 

I*>* = E<pIp>- (83) 

p=0 

This basis set is infinite due to the non-compactness of u(l, 1). The basis states 
are generated by a raising operator Z + 

\p) = ^rl"). (w) 

where is the lowest basis state. The basis states are orthogonal with norm 

As n — > 0, (p\p) = 5 0p . Bocquet's Z is identical to the Z' in this paper. 

We use the same general form for the ground state of our transfer Hamiltonian 
but summed over an additional variable 

oo 

l*>*= E C PP '\P,P')- (86) 

p,p'=0 

The non-primed terms are related to the advanced Green's function and the 
primed terms are related to the retarded Green's function. The primed term 
should act in the same way as Bocquet's variables and the non-primed should 
act like complex conjugates. 

We construct some raising and lowering operators, Z± = Z y =p iZ x and Z' ± = 
Z' y TiZ' x (as was done in section 5). The basis states, \p,p f ), are generated by 
the action of Z' + 

\p, P ') = { -^\p,n'), (87) 

and, since the non-primed terms act like complex conjugates of the primed 
terms, the action of Z_ 

(Z ) p 

\p,p')= { —^\{l,P>). (88) 
p\ 

The state \ is the lowest vector in this representation. Alternatively, we 

could have Z + generating the basis states since [ZS)^ = Z + 

(p,p'\ = (n, P '\^. (89) 
p] 
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The vector Z' operates on the \p,p') basis state in the following way 

Z' z \p,p'} = (p' + f) \p,p') 
Z' + \p,p') = (p' + l)\p,p' + l) 

Z'_\p,p') = (p' + n'-l)\p,p'-l), (90) 

and for Z operating on (p,p'\ 

{p,p'\Z z = (p+ f) {p,p'\ 
{p,p'\Z+ = (p + l)(p+ l,p'\ 

{p,p'\Z_ = (p + n-l){p-l,p'\. (91) 



The scalar products of these states can be calculated by using, for example, 

(z'_\p, P ')y = (p, P '\z' + 



{p,p' - l\(p' + n' - 1). 



(92) 



The norm is 



(p,p'\p,p') = (p,p 



-Az_y(z' + r' 



p\p'\ 

T{p + n)T{p' + ri) 



T(n)T(n')T(p+ l)T(p> + 1) 
T(p + n)T(p' + ri) 



(n,tf\n,ct!) 



T(n)T(n f )T(p+l)T(p' + 1)' 
Note that Z and Z' commute. In the limit n, n' — > 

{p,p'\p,p} = 5 p5 p>. 



(93) 



(94) 



Similarly, it can also be shown that the scalar product of \p,p') and \q, q') with 
p ^ q oi p' ^ q' vanishes, assuming {p,p'\fl,Q') = if p ^ Q and p' ^ Q'. 
Hence, the basis states are orthogonal. 

As discussed previously, because the transfer Hamiltonian isn't Hermitian the 
left wave function is not the Hermitian conjugate of the right wave function. 
The Hamiltonian and its Hermitian conjugate are related by the transform 
R = ex.p(iir(Z z + Z' z )) [7]. The left and right ground state wave functions are 
related by this same transform 



\*)l = R\*)r. 

Using equation (93) and the orthogonality of the basis states 
mm\ 1 y ( lY+P ' + ^ nP + n)T(p' + n') 



(95) 



(96) 
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As n, n' — > this reduces to 

L(*m R = ICool 2 (97) 

so we can make our normalization condition | Coo 1 2 — 1, m agreement with 
equation (66). 

Using equations (90), (91) and (96) 



( n'\ T{p + n)T{p' + n') ^ 2 
Thus, by substitution into equation (67), 

(G A (E,x,x)^G R (E,x,xy) = ^mj) ? ( - 1)P+P ' 

Xp k-j-i p/j -i\ Cppi \2 (9Q) 

If we can find the ground state of the transfer Hamiltonian we can calculate 
all the C PP ' and hence we have calculated the Green's function correlator. We 
will calculate ( pp r for a special case of the Dirac equation corresponding to a 
single channel quantum wire in section 7.1. 



6.2 Super symmetry 



The bosonic parts of the superspins, J and J', are defined in exactly the same 
way as the Z and Z' used in the replica trick calculation when n — n' — 1. So, 
J and J' satisfy all the same equations as Z and Z' when we set n — n' — 1 in 
all the replica trick calculations. Like the replica trick we have a bosonic wave 
function \p,p') but here we will use the notation \p)\p')' . This bosonic wave 
function is not affected by the fermionic operators S and S'. From equations 
(90) and (91) 

JM={p + \)\p)\ ( P \J Z = ( P + 1) (pi 
J+\pY = (p + i)\p+ i)', (p\J+ = (p + i)<p + 1|, 

j'_\ P y = p\p-iy, (p\j- = p(p-i\, (ioo) 

where J± = J y ^f iJ x and J' ± = Jy =)= i J' x . From equation (93) 

(p\p) = '(p\pY = i. (ioi) 
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The fermionic operators, S and S', are like spin operators. There are two spin 
up states and two spin down states and the fermionic operators operate on 
them as follows 



s;ir>' = 


w m 


(T \S. 


= 1(11, 


S' + \ T>' = 


o, 


(T \S + 


= 0, 


S'_\ T>' = 


-1 1)', 


(T \S- 


= -ai 


s' z \l)' = 


— 2I J-) ' 


U \S, 


= -ka 


s' + \l)' = 


1 T>', 


(i\s + 


= <H 


SL\i)' = 


0, 


u \s- 


= 0, 



(102) 

where S± = S y =fiS x and S'± = S' y ^iS' x . Note that the operation of S- and S'_ 
on the spin up has the opposite sign from usual. This is due to the Grassman 
variables which make up the spin operators. For the same reason it can be 
shown that (S z ) 3 = (S' z ) 3 = 0. 

Balents and Fisher [6] show that the basis state for the non-primed bosons 
and fermions is 

,'b>ll) + |p-i>IT>, p>o, 

\p)r = , v , 4 (103) 
|0)||), p = 0, 

and similarly for \p)' R . The total basis state combines primed and non-primed 
parts and is written as \p)r\p)' r - The left basis state for \p) R is 



\p)l = (~l) J ' 



R 

(-l)p(|p)|i) + |p-l)|t)), p>0, 
|0)||), p = 0, 



(104) 



and similarly for \p)' L . 



As discussed in Bocquet [7] the norm of the lowest basis state, \0)r, can be 
taken to be unity 

L <0|0) fl =<0|0)UU) = l. (105) 

Since the |p)'s are orthonormal the norm for the down spins must be unity. 
The up spin can be defined in terms of the down spin: 

(T I = (I \S+. (106) 
This allows us to write the norm of the up spin as 

<T I T> = (I \s+S-\ t), (107) 
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which can be evaluated using the spin rules of equation (102). We can calculate 
the norms of the primed spins in the same way. The norms for all the spin 

states are 

-'<T|T>'= , UU> , = -<TIT> = aU> = i- (108) 



The total left and right wave functions are created by summing over all the 
basis states 

oo 

I*>l = E C pp >\p)l\p')' l , 

p,p'=0 

oo 

l*>*= E C pp >\p)r\p'Y r ,. (109) 

p,p'=Q 

Using the rules given above and the normalization condition on the wave 
function it can be shown that 

L (tt|tt) fl = iCool 2 = 1. (HO) 



Combining equations (100) and (102) we obtain the following equations for J 
and J' 



J'zWr = p\p)'r, l(p\Jz = l(p\p, 

J' + \p)'r = (P + l)b + 1>h. l( P \J + = l{p + l|(p + 1), 

J-\p)' R = (P - i)b - 1)r, l{p\J- = l( P - l|(p - i). (in) 

These are exactly the same as Z and Z' for the replica trick for n = n' = 0. 
To calculate (79) we require 

L (*l = m{j z - s x y-*{j' z - s>y\*) R . (ii2) 

This expression can be expanded using equations (102) and (111) and by using 
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the form of the ground state in equation (109) 

Jt j J'jm R = (J z - s.y-'w - s' z y\*) R 

= [Jt j - {k- 3 )jt 3 - l s z + Uk-m-j - i)j z k - j ~ 2 si 

x (J'i - jJ't X S' z + - l)Jt 2 S> 2 ) \9) R 

oo 

E (p k ~ j - (* - jtf-^s, + \{k- j){k - j - i) P k - j - 2 sl) 

p,p'=0 

x (p'l-w'^S: + - l)j*-*Sf) Cpp>\p)r\p')' r 

oo 

= E l(p' j + W~ x + - l )p' j ~ 2 ) W)'\l)' 

p,p'=0 

+ ( P ,j - w 3 ~ l + yu - i)p°" 2 ) \p' - i)i t)'" 

x ~(p h - j + \{k - 3)p k -'- X + \{k-j){k-j - l)p k -i- 2 ) \p}\ I) 

+ (p k ~ j - \{k-j)p k ^~ l + l(k-m-j - iy^'- 2 ) 

xb-i)l T)]C PP ', (113) 

where, in order to simplify the above formula, we have set | — 1)r = | — 1)'^ = 0. 
Using equations (108) and (104) and the orthogonality of the basis states 

oo 

L <*|j;-^i*>fl= e {-iy +p '\CpA 2 {k-i)jp k -i-^ 
p,p'=i 

+ E(-l)1C P o| 2 (A: - J)p k ~^ (p'i + W 3 ' 1 + - 1)^'~ 2 



P =i 



p'=0 



+ E(-l) P 'Kopf jp' j ~ l (p k - J + \{k - j)p k -i~ l 
P '=i 

+l(k - j)(k - J - l)p k -^ 



p=0 

+ (p k - J + \{k - 3)p k -i- 1 + \ik — j)(k — j - l)p k ~^ 2 ) p=Q 

x ( P ' j + \tf'- x + bU - i)^'~ 2 ) p , =0 • (H4) 

The terms arising fromp = will vanish for k—j > 2. The terms due to p' = 
will vanish for j > 2. The above equation may be substituted into equation 
(79). When k—j > 2 and j > 2 this equation will resemble equation (99) 
but with the Z z and Z' z replaced with J z and J' z . For k — j < 2 or j < 2 the 
supersymmetric form will have some additional constant terms. However, as 
will be shown in the following section, the density of states calculated from 
supersymmetry does reduce to the density of states calculated from the replica 
trick. As for replica trick we need to solve for ( pp > in order to calculate the 
Green's function correlator. 
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7 Evaluation of the moments of the local density of states 

We shall now evaluate the moments of the local density of states for a specific 
case corresponding to a disordered quantum wire with a single channel. We 
set = m = v = D v = and = D m = D. 

7. 1 Solution of the ground state wave function of the transfer Hamiltonian 

The replica trick transfer Hamiltonian (57) with <f) = m = v = D v = and 
= D m = D is 

H = 2iE(Z z - Z' z ) + 2 V (Z Z + Z' z ) 

+ 2D {{Z + - Z' + ){Z_ - Z'_) + (Z_ - Z'_)(Z + - Z' + )) . (115) 

The supersymmetric transfer Hamiltonian (76) has a similar form but with the 
vectors Z and Z' replaced with the superspins J and J 1 . Also, as can be seen 
by comparing equations (90) and (91) with equation (111) the operation of Z 
on the replica ground state with n = n' = is identical to the operation of the 
J on the supersymmetric ground state. Hence, the ground state eigenvalue 
equations must be identical when we take n and n' to zeroth order in the 
replica ground state. We shall calculate the zeroth order replica ground state 
here and use the same ground state in the supersymmetry calculations. 

If substitute the transfer Hamiltonian into 

(p,p'\H\*) = 0, 

using, for example, 

(p,p'\V) =(pp> 
(p,p'\Z z \V) =p( pp/ 

{p,p'\Z + Z^)=(p + l)p( pp , 
(P,P'\Z+Z-\V) =PP'(p+ip'+i, 

we obtain 

iE(p - p')( pp , +r](p + p')( pp/ 

+ 2D(p 2 ( ppl + p' 2 C Ppl - pp'(p-ipi~i - PpXp+ip'+i) = 0. (118) 

This equation implies that Q p , = ( p ' p which implies that ( pp is real. If we 
assume that ( pp r = d~ pp '( pp then 

VPCpp + Dp 2 (2( pp - Cp-ip-i - Cp+ip+i) = 0. (119) 
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(116) 



(117) 



Note that this equation is the same as equation (47) in Reference [19] which 
was obtained using the Berezinskii technique. We define a new variable x = 
rjp/D. As p increases to p + 1 x only increases by rj/D. It is reasonable to 
assume that r] <C D so we can take a; to be a continuous variable. If ( pp = ((x) 
we obtain the following differential equation of ( in terms of x 



d 2 ((*) _ I c(x) 
dx 2 x 



0, 



(120) 



which has the general solution 



C(x) = x 1 ' 2 la'K^x 1 ' 2 ) + bh{2x^ 2 ) 



(121) 



where a' and b' are constants. In terms of p 



5PP 



p 



1/2 



aKl[2 m +u j,m 



(122) 



where a and b are constants. 

Boundary conditions require that ( remains finite for all p so we must have 
6 = 0. We can find the other constant of integration, a, by the normalization 
condition on the ground state. The normalizing condition is | Coo 1 2 = 1- We 
also know that Coo must be real so 



: 00 = l = a\imp^K 1 (2^ 

p^o V V D 



(123) 



which gives 



a = 2./£ 



(124) 



Hence 



(125) 



and 



Cpp' — dpp'^^I^Ki ^2 




(126) 



28 



7.2 Replica trick 



By substituting equation (126) into (99) 

{2i)\-iy (D\ k - 1 k T(k)* 



r(k-j)r(j) \ V J fc-ir(2fc)- 

We have used the fact that in the limit i] ^ the sum over p can be con- 
verted into an integral by defining the variable x = 2^Ji]p/D which may be 
approximated to a continuous variable since r\ <C D 



p=0 



'dV' 1 k r(k) 



k- 1 r(2fc) ' 



(128) 



and [34] 



r^^^'i^ifcl)' (129) 

if A > 1. Note that as rj — > the sum in equation (127) is dominated by 
the terms with large p. Hence, the non-compactness of u(l,l) is playing an 
essential role here. 

By substituting (127) into equation (25) we obtain the moments of the density 
of states 

(Pl) = ^UJ (2FTT) (130) 



by using the relation (calculated by Maple 5.1) 

^ (k-j)j (k-r 
f ((k-j)\j\y r{k) 



,,f (k-j)j _ (k - i)m - 1) 

ft -z^ rm -MAM - wtaz • ^ L6L > 



When k — 1 we find that (pi) = 1/ir. This equation for (pj) agrees with 
equation (69) in Al'tshuler and Prigodin [19] and, more recently [21], where 
Si = 8i]r = i]/D. The two results differ by a factor of l/n k because what 
Al'tshuler and Prigodin define as (p k ) we define as (p k ) / (p) k , and similarly for 
all the density of states. 
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7.3 Supersymmetry 



To calculate the moments of the density of states from supersymmetry we 
substitute equation (126) into equation (114). The terms from p = and 
p' — will be negligible for small r\ for all cases except k — 1. From equation 
(114) and for k > 1 we obtain the same result as for the replica trick, shown 
in equation (127). The moments of the density of states are given in equation 
(130). 

The calculations for k — 1 are straight forward. It is easy to show 

l (*\Jz\*)r= L(nm) R = i (132) 

which gives (p\) = - which agrees with equation (130). 



8 Moments of the full density of states 



We showed at the end of section 5 how to calculate the moments of the full 
density of states, (pi), described by equation (82). We require these moments 
rather than the moments of the density of states which removes the rapidly 
varying terms, (pi), described by equation (80). By comparing these two equa- 
tions we see that to obtain the moments of the full density of states instead 
of J'J we need the term 



z 



3 m -| 2ik F (l+m-j) 

hU~ 2 ' + " ' u -my.( m -,y.,<: (133) 

A similar expression may be written for the additional non-primed terms. We 
have ignored that fact that the components of 3' do not commute since, as will 
be shown, the terms arising from interchanging the components are negligible. 
Instead of J?~ j J'J we use 

3 k—j m r 

( 1 ) k-m-r+l+q jij—m jll j/m—l jk—j—r jq jr-q 

m=0 r=0 1=0 q=0 

x tlZ J ±l . (134) 

(j — m)\(m — l)\l\(k — j — r)\(r — q)\q\ 

Now we can neglect those terms which contain rapidly varying exponentials. 
We retain those terms for which q = k — I — m — r. However q must be in the 
range [0, r]. This places a restriction on I 

max(0, k — m — 2r) < I < (to, k — m — r). (135) 
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The above sum is now 

j k—j mm(m,k—m—r) 

•y •y ^1^2(k-m-r) j'j-rn j/l jfm-l jk-j-r jk-m-r-l j2r+m+l-k 

m=0 r=0 i=max(0,fc-m-2r) 

J'K*-J')! 



(j — m)!(m — /)!/!(& — j — r)!(2r + m + / — — m — r — l)\ 



. (136) 



We need to calculate something of the form l(^\J+J-J'zJ+J-JI\^)r- To get 
insight into how this is done we first evaluate 

oo 

L <*| J$J^J?\*) R = £ \C PP \\p- 1) • • • ( P -b + l)(p-b + l) . . . (p-b + a- 1) 

x [p(p - b + a)(p c + \cpc~ 1 + |(c - l)(c - 2)p c " 2 ) 

- (p - 6) V - ic^" 1 + |(c - 1)( C - 2)p- 2 )], (137) 

where we have used ( pp > = S pp '( pp . As rj — > the sum is dominated by terms 
with large p (since C PP is proportional to 77) and so we retain just those terms 
with the highest power of p which is p a + b + c - 1 . This is why in equation (133) we 
can swap the components of J' around, because the extra terms are a smaller 
power of p than p a+h + c ~ l which become negligible as rj — > 0. The simplified 
result is 

00 

l(*\J!?J?J?\*)r = E + b + c)p a+b+c -\ (138) 

p=0 

When we include both the primed and non-primed terms we obtain 

00 

L <* I J" J* J? J* Jl jf\*) R = J2 \C PP \ 2 (a + b + c)(d + e + f)p a + b + c +d+ e + f-2_ 

(139) 



Substituting these identities into equation (82) gives 

" p=0 j^J-y 1 ^ JJ- m =0r=0 

mhi(m,k—m—r) 

x J2 (|) 2 ( fc - m - r ) 

i=max(0,fc— m— 2r) 

x [(j - m)\(m - l)\l\(k - j - r)!(2r + m + I - k)\(k - m - r - Z)!] _1 . (140) 

The sum over p can be evaluated as in the previous section using equation 
(126) and (128). 
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The sum over j, m, r and I can be rewritten as [23] 

^ = E4^4yE EE E Mi) 2(fc - m - r) 

j=0 J ' v J ) ■ m=0 r=0 Z=0 l'=k-m-2r 

x [(j - m)!(m - l)U\(k-j - r)\(2r + m + I' - k)\(k -m-r- Z')!] -1 . 

(141) 

If we write the kronecker delta as 

Sw = 7T / (142) 
2n jo 

and rewrite the sum over /' as a sum over q = I' — k — m + 2r it can be shown 
that 

s= (fc-l)r(2fc-l) 1 j -a + 1 + e « )fcj (143) 



fcr(A;) 4 2tt 7o 

by making use of equation (131). The integral can be evaluated after a binomial 
expansion. After integration the resulting sum may be simplified using Maple 
5.1 so that 

(fc-i)(2fe-i)r(2fc-i) 2 

2 fc -W(A;) 6 • 1 J 



The moments of the total density of states is thus 

This agrees with equation (68) in Al'tshuler and Prigodin [19] although we 
have an additional factor of l/7r fe , as discussed previously, due to different 
scaling of the density of states. Using equation (23) we can use the moments 
of the density of states in equation (145) to calculate the fcth moment of the 
participation ratio 

P (k) = {2D) k - 1 T{k), (146) 
where we have used (p) = lim 7? _ > o(pi) = 1/tt. 



9 Probability distribution function 



The moments of the participation ratio may be used to find moments of density 
of states, (p k f), which are derived from some arbitrary weight / as defined in 
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equation (24). The probability distribution function of any local density of 
states, W(pf), may be obtained from these moments of the density of states 
by the following [19] 



/ W(p) exp(—pp)dp = exp 
Jo 



The cumulants are defined by 



£ 

k=i 



-p) k (p k f ) 

kKPf) k 



(147) 



( P ))c = ( P )) - ( P k f U 

Solving equation (147) for W requires an inverse Laplace transform. 



(148) 



Following Al'tshuler and Prigodin [19] we now summarise how one can con- 
struct a density of states which is related to the NMR line shape. In a metal 
the NMR line undergoes a frequency shift, known as the Knight shift, which 
is proportional to the local density of states derived using (5) with weight 



f2(E) 



n F (E + gp e H) - n F (E - gp e H) 
2gp e H 



(149) 



This weight is the derivative of the Fermi distribution function, n F , about the 
Zeeman splitting gp e H. The Fermi distribution function is 



n F {E) 



exp 



E — e F 
T 



-i 



(150) 



where e F is the Fermi energy. The density of states derived from this weight 
within the region T ^> gp e H will be defined as p2(T, e F ). In this high temper- 
ature limit 



/ 2 (£) 



1 

4T 



cosh 



-2 



E - e F 
2T 



(151) 



The local density of states, p 2 , fluctuates throughout the sample. Hence, the 
Knight shift also fluctuates. The net result of these various Knight shifts is a 
broadening of the NMR line shape. If the relaxation rate is of the sample is less 
than the Knight shift the NMR line shape is fully determined by the super- 
position of all the local Knight shifts. This is proportional to the distribution 
of the density of states, p 2 . 

By substituting equation (146) and (151) into equation (24) we find the mo- 
ments of this new density of states are 



(P2(T,e F )-)= — 



277 



1 r{k) 2 r( 



7rr(A; + i)' 



(152) 
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where we have used [34] 



dy cosh 2k y = 2 



T{k)T 



T(k + 



(153) 



The sum in equation (147) is dominated by the large k terms. If k is large we 
can assume that the /cth cumulant is approximately equal to the kill moment 



1.5 



1 - 



0.5 r- 








Fig. 1. The probability distribution function, equation (155), for a number of dif- 
ferent values of s = 2T/irD. For large s (high temperature and low disorder) the 
distribution is symmetric so that the typical value of the density of states is the 
same as the average value. As s decreases (low temperature and high disorder) the 
distribution becomes skewed and there is a marked difference between the typical 
and the average value of the density of states. This probability distribution can be 
related to the NMR line shape. 
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of the density of states 



<P2> C = <P2>- (154) 

After substituting the moments shown in equation (152) into equation (147) 
and taking the inverse Laplace transform Al'tshuler and Prigodin [19] obtained 



s f°° 

W{p2) = — / <ir sinh2:r sin(7rs:r) exp 

Z7T J-oo 



P2 , 2 2 71 

■ cosh x + x — - 



(P2> 4 



(155) 



where s = IT j ttD in our notation but s = 16Tt/it in Al'tshuler and Prigodin's 
notation. The relationship between the disorder strength, D, and the scatter- 
ing time, r, is 8r = D^ 1 . 

The probability distribution function at large temperature of the density of 
states, p2, is plotted in figure 1 for a number of values of s. Note that for 
s < 1 the distribution is quite asymmetric and the typical (most probable) 
value is significantly less than the average value. We have set various constants 
to unity, such as the Boltzmann constant, ks, the Fermi velocity, vp, and h. 
If all these constants are included we find that near the Fermi energy s = 
2kBThvF/nD or, in terms of the mean free path, / = vft = {Tivf) 2 /SD, we 
have s = l^ksTl/^TiVF = 6.67 x 10 11 TI/vf- If we have a sample with a mean 
free path of 1 pm and a Fermi velocity of vf — 1 x 10 6 ms _1 then some skewness 
in the NMR line shape should be observable when the temperature is 1.5 K or 
less (s < 1). 



10 Conclusion 



We have presented a general formalism to calculate exactly the disorder av- 
erage of a product of Green's functions, (G A (E,x,x) k ~ : 'G R (E,x,x) : '), for the 
random Dirac equation using both the replica trick and supersymmetry. This 
extends the formalism developed recently to calculate the disorder average of 
one Green's function using the replica trick [7] and supersymmetry [6,7]. The 
problem is reduced to that of finding the ground state of a zero-dimensional 
Hamiltonian involving a pair of "spins" that are elements of u(l, 1). The mo- 
ments of the different density of states can be written in terms of the expec- 
tation values of various "spin" operators in this ground state. We explicitly 
found this ground state wave function and the associated expectation values 
for the case of the Dirac equation corresponding to a disordered quantum 
wire with a single channel. The non-compactness of u(l, 1) played an essential 
role in the results. We showed how the moments of the local density of states 
changed when it was averaged over atomic length scales. Our results for the 
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moments agree with those obtained previously by Al'tshuler and Prigodin [19] 
using the Berezinskii diagram method. In our view, our derivation is more 
transparent from a field-theoretic point of view. Future work should focus on 
finding the ground state of the transfer Hamiltonian corresponding to more 
general random Dirac equations, particularly those associated with quantum 
phase transitions in random spin chains. 
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A APPENDIX 



In this appendix we derive the prefactor in the replica trick expansion of the 
kth power of the Green's function shown in equation (33). Consider an N 
dimensional integral of the form 



IN 



J d N xF(x 2 l ) (A.l) 



where F( function which depends only on the length r = YliLitf- By 

making a change of variables into polar coordinates: 

(x!,x 2 , ... , x N ) -> (r, (f), 0i, . . . , 0Ar_ 2 ) (A. 2) 

we can evaluate the integrals over the 0's and the <fi (since F(xf) = F{r) has 
no dependence on the 0's or the 0) and we are left with [35] 

n N/2 roo N _ 2 

irr ^ F(r) - (A ' 3) 



For example, if s and s* are one dimensional variables such that s = x\ + ix 2 
and A is a constant 



1 



1 



* „—s* As 



ds*ds(s*s) K e 



ds*dse~ s * As 



i -l 



(A.4) 



r k e Ar where r 



s s = 



then the integral has dimension N = 2 and F(r) 
x\ + x\. Using equation (A.3) it can be shown that in terms of r the integral 
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IS 



1 



1 



/ drr k e~ Ar 




/ dre~ Ar 


Jo 




Jo 



(A.5) 



We now consider a replicated version of equation (A. 4). If this equation is 
replicated n times 



1 X 



„ n / n \ k / n \ 

y n rfs /* rfs « s ^J ex p z s *^J 

/n /it 
Y[ dsldsi exp I - ^ s*As t 
i i \ 7 i 



(A.6) 



where X is a factor which must be calculated. We can write sj = X21-1 +1x21 so 
now we are dealing with a N = In dimensional integral and r = J2?=i s i s i — 
Ya=iXi- The integral becomes 



J_ X 7r"/r(n) 
A* ~ fcl n n /T(n) 



n—l„k„—Ar 



drr n re 



1 

Jo 



drr n-l r -Ar 



(A.7) 



The integrals over the r variable may be evaluated to show that 

J_ X 1 T(n + A;) 



so 



k\F(n) 



T(n + k) 

and hence the replicated form of equation (A. 4) is 



(A.8) 



(A.9) 



1 



T(n) 



A k r(n + fc) 



/n / n \ ^ / n > 

/n / n 
Y[ dsldsi exp - ^ s z *As; 
7 1 \ 7 1 



. (A.10) 



Note that when n = 1 we get back to equation (A. 4). If k — 1 then the factor 
before the integral is 1/n which agrees with Bocquet [7]. 
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